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ABSTRACT 


Temperature studies vere conducted during both simulated and human 
occupancy tests on the NRDL 100 man shelter. The temperature of the 
shelter remained comfortable during l4 days of human occupancy. The 
ventilating air movement through the shelter designed for l6 cfm per 
occupant was adequate to reduce the carbon monoxide concentration from 
heavy smoking to below safe tolerance limits and it removed 69 # of the 
heat generated by the occupants with only 31 $ of the heat being dissi¬ 
pated through the walls. The shelter temperature varied between 7^ and 
90°P with an outside air temperature of 50 to 94°F during the simulated 
occupancy test and between 70 and 82°7 during the human occupancy tests 
with an outside air temperature of 36 to 67 °P. Poor distribution of 
the inlet air resulted in uncomfortably cool areas In the front of the 
shelter during the human occupancy test. 

A preliminary analysis was conducted on a thermal analyzer to 
determine its potential as a tool for predicting shelter temperatures 
under various climatic conditions. This analysis showed only about 
10 $ deviation from the measured temperature reported here. 
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SUMMARY 


The Problem 

To determine the adequacy of the NRDL 100 man shelter design to 
maintain comfortable living conditions for at least l4 days of human 
occupancy. 

Findings 


Under the conditions of the test the temperature and humidity of 
the shelter remained within comfortable limits for 14 days of human 
occupancy test except for cool drafts at the front of the shelter due 
to poor distribution of the inlet air. The average heat output by the 
occupants was 485 Btu/hr per man. The ventilating air movement through 
the shelter of 16 cfm per occupant was adequate to reduce the carbon 
monoxide concentration due to heavy smoking to below safe tolerance 
limits. It removed 69 of the heat generated by the occupants with 
only 31 of> the heat being dissipated through the walls. The effective 
temperature in the shelter varied between 74 and 80° during the simu¬ 
lated occupancy run when the ambient air temperature varied between 50 
and 94°F. During the human occupancy test the effective temperature 
was 70 to 82° with the outside air temperature varying between 36 and 
67°F. 


A preliminary study on a thermal-electrical analyzer at Oklahoma 
State University showed approximately 10 $ variation from the measured 
conditions. This indicated that reliable studies can be made in the 
laboratory for predicting the environmental conditions in this shelter 
located In various climates. 
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CHAPTER 1 


INTRODUCTION 


The protection of the nation from the hazardous radioactive fall¬ 
out of a nuclear war has been given much consideration in recent years. 
The most common considerations are the protection of individuals in 
family shelters and the protection of large groups in community-built 
shelters. 

A prototype 100-man shelter was constructed by NRDL in 1959 at 
Camp Parks approximately 40 miles east of San Francisco. A variety of 
experiments have been performed in this shelter to establish adequacy 
of the ventilation in maintaining comfortable temperature and humidity, 
to study the psychological effects of confinement in a group shelter, 
to study the need for closure during mass fires, and fallout ingress 
experiments to determine the need for filters in the ventilation system. 
The results of these studies are abstracted in reference (l). 

In the design of a community shelter it 1 b essential that the cost 
be kept to a minimum, whether these shelters are to be financed by the 
government or by private enterprise. The problem of building a shelter 
at lowest possible cost automatically requires that both the floor area 
and the gross volume per man be reduced to a minimum. To maintain 
physical comfort in a shelter under these conditions, sufficient venti¬ 
lating air must be provided to carry off the heat and moisture generated 
by the occupants. These factors were weighed against cost in the design 
of the U. S. Naval Radiological Defense laboratory 100 man shelter and a 
blower was selected which was expected to maintain a comfortable envir¬ 
onment at minimum cost.2 

The methods available for theoretical analysis of the heat-transfer 
conditions in the shelter were limited so it became necessary to con¬ 
duct temperature studies during actual human occupancy tests. As a 
precautionary measure prior to conducting these tests a study was made 
with 100 devices which simulated human occupants. 
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CHAPTER 2 


EXPERIMENTAL CONDITIONS 


2.1 SHELTER CONSTRUCTION 

Figure 1 is an artist's sketch showing the basic design features. 

It is a galvanized multi-plate steel arch type structure (Standard Navy 
Ammunition Hut) with a floor area of 25 ft x 48 ft. It is burled with 
the top flush with the ground level and then covered with a 3 ft mound 
of soil. It has a capacity of 100 people with a floor space of 12 sq. 
ft per man and a gross volume of 117 cu ft per man. It has 4 tiers of 
bunks on each side with a sleeping capacity of 96 . These bunks can be 
rolled up to provide more space in the daytime. Air intake openings 
are located in the entry-way at the top of the stairs as shown in Fig. 

2. Fresh air is drawn in with a blower located outside the inner door. 

The air comes into the shelter space through a duet over the inner door 
and is exhausted out a vent located at the rear. The air intake can be 
closed and the exhaust vent dropped down to give an air tight closure. 
Tanks are provided for water and gasoline storage, the gasoline being 
used to operate an emergency power unit. 

The soil surrounding the shelter was of two types as shown in Fig. 

3 . Construction specifications called for backfilling up to a level 8 ft 
from the floor with select fill compacted to a specific density. The con¬ 
tractor chose to use pea gravel for this select fill. The original 

soil removed during excavation was then put on top of this select fill. 


2.2 INSTRUMENTATION 

2.2.1 Temperature Measuring Stations 

Upon completion of the shelter construction, holes were drilled 
through the walls and thermocouple probes were inserted at the loca¬ 
tions shown in Fig. 4. 

A series of probes were located on the center line at 4 ft, 8 ft, 
and 11 ft from the floor, at the top of the arch and one in the floor at the 
center. An additional set was located midway between the front and the 


2 



AIR INTAKE 
VENT MUSHROOM 


EXHAUST 

VENT 


EXTERIOR DOOR 



Fig. 1 Cut-avay Diagram of Shelter 
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Fig. 2 Shelter Entrance Shewing Blast Bulkhead, Exterior 
Door, Mushroom Shaped and Intake Vent Covers 


4 






SURFACE AREA OF SHELTER IN CONTACT 
WITH THE VARIOUS SOILS 

SOIL CONTACT AREA 

SELECT FILL (SF) 1160 SQ FT 

BACK FILL (BF) 1175 « « 

FLOOR (F) 1200 " " 


Fig. 3 Soil Conditions Surrounding the ITRDL Shelter 
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B45T 



F- FRONT SECTION L-LEFT SIDE 

C-CENTER LINE T-TOP 

B- REAR SECTION R-RIGHT SIDE 

NUMBERS DESIGNATE THE HEIGHT 
ABOVE THE FLOOR. FOR EXAMPLE- 
LF-B-LEFT FRONT 8FT FROM 
FLOOR. RC-11-RIGHT CENTER 11 FT 
FROM FLOOR. 

F45T 1 THERMOCOUPLE LOCATED 
&B45TJ IN CORNER AT 45° ANGLE. 


Fig. L Temperature Measuring Stations 
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center line and another one midway between the rear and the centerline at the 
8 ft level. A single probe was located in the back wall at the 8 ft 
level. The thermocouples at each station consisted of five couples, 
one on the inside of the metal surface, one at the outside surface and 
one each at 1 ft, 2 ft and 3 ft into the surrounding fill. 

The external surface and the thermocouples at 1 ft, 2 ft and 3 ft 
from the shell of the shelter were positioned on thin-wall, half-inch 
diameter phenolic resin tubes that were pushed through holes in the 
shelter wall (Fig. 5)* The thermocouple leads were secured to the out¬ 
side of the plastic tubes at proper locations. After positioning the 
four couples on each tube, one end was sealed and the tube then filled 
with sand. A half inch pipe with a loosely fitting steel point in the 
end was driven into the soil at the desired locations and the probe 
slipped into place. The pipe was then pulled out over the lead wires, 
leaving the steel point and the probe in the soil. 

The internal surface temperature thermocouples were attached to 
the wall with masking tape. Wet and dry bulb temperatures were taken 
a\. the entrance, center of shelter and at the exhaust. The thermo¬ 
couple leads were all taken out through an opening in the top of the 
shelter near the exhaust vent and connected to continuous strip chart 
recorders. Continuous temperature recordings were made with an auto¬ 
matic printout of each temperature every 8 minutes. 

During the simulated occupancy tests, wet bulb thermocouples were 
located at the entrance to the shelter, inlet over the inner door, cen¬ 
ter of the room and at the exhaust. The evaporation at these thermo¬ 
couples was accomplished hy a small electric fan. 

2.2.2 Human Simulators 


The simulators were designed to simulate the latent and sensible* 
heat output of a human. These human simulators consisted of a nichrome 
wire heating coil mounted in the bottom of a five gallon paint pail. 

The pails were covered with cotton flannel to serve as a wick. The 
moisture output of the simulators was provided by dripping water onto 
the top of the pail. The water output was regulated by a valve at each 
simulator and the total flow of water into the shelter was controlled 
in the instrument shack. The amount of water to be added was varied 
with the temperature at the center of the room. 

■^Sensible heat added to the shelter is a direct addition of heat to the 
enclosure by any cneor all of the mechanisms of conduction, convection, 
and radiation. A gain of latent heat is considered to occur when there 
is an addition of water vapor to the air of the enclosure. 
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SHELTER 



Fig. 5 A. Pipe Driven Into Soil Prior to Inserting Thermocouple Prohe 
B. Thermocouple Prohe in Location in Soil After Removing Pipe 











Thirty-six simulators were located on each side in the approximate 
position of the center bunk as shown in Pig. 6. , 

The remaining 28 simulators were located in the front section at 
the tables as shown in Fig. 7* 

2.2.3 Carbon Monoxide Measuring Equipment 

The carbon monoxide measuring instrument was provided by the 
U. S. Forest Service. This equipment, developed by the Forest Service 
Pacific Southwest Forest and Range Experiment Station, consisted of 
surplus detectors which had previously been used as warning device in 
aircraft. These detectors split the flow of a gas sample through two 
cells connected by a thermopile. One cell contains activated hopcalite* 
which catalyzes the oxidation of carbon monoxide. The resulting temper* 
ature difference between the two cells is a function of the carbon mon¬ 
oxide concentration in the air being sampled. This temperature differ¬ 
ence results in an emf being generated in a thermopile which is recorded 
on a strip chart recorder. The unit is calibrated to give approximately 
100 millivolts for a concentration of one percent (by volume). 


*Hopcalite - a mixture of oxides of copper, cobalt, manganese and silver, 
used in gas masks as a catalyzer converting carbon monoxide to carbon 
dioxide. 
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CHAPTER 3 


TEST PROCEDTTRES 


3.1 TEMPERATURE TEST'S 

A total of five tests were conducted as shown in Table 1. During 
the simulated runs No. 1 and 2, the water output of the simulators was 
regulated by a valve at each station and the total flow of water into 
the shelter was measured by a flcrwrator in the supply line. The amount 
of water added was varied with the center room temperature in accord¬ 
ance with the values given in Pig. 8. The total heat into the shelter 
by the simulators was maintained at 4?0 Btu/hr per man, which is an 
estimate of the heat that would be generated by a shelter occupant. 

This value is based on data from the Heating and Ventilating Guide of 
the American Society of Heating, Refrigeration and Air-Conditioning 
Engineers which give 3^0 Btu/hr as the heat generated by a human seated- 
at rest, 430 standing at ease, and 760 , walking at 2 mph. 

Preliminary heat balance calculations on runs No. 1 and. No. 3 indi¬ 
cated inaccuracies in the heat transfer coefficient for soil, because 
of the diurnal temperature variations in the inlet air. Therefore runs 
4 and 5 were conducted with air vents closed. These tests were continued 
for a sufficient length of time to establish a constant rate of heat loss 
to the soil. 


3.2 METHOD OP PREDICTING SHELTER THERMAL CONDITIONS 

Heat transfer in this shelter is a complex process and predicting 
the shelter environment under various climatic conditions is difficult. 
Therefore arrangements were made with the Mechanical Engineering Depart¬ 
ment of Oklahoma State University to determine the feasibility of using 
a system which they have developed and called a thermal analyzer to 
determine shelter thermal behavior under various specified conditions. 
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LBS WATER/HR ADDED TO SHELTER / MAN 


Fig. 0 Relationship of Water Emitted, by 100 
Persons to Temperature of Their Environment 
(based on data from the ASUVE Guide) 
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The thermal analyzer applies the principles of thermal-electrical 
analogy which is based on the similarity of the equations describing 
heat and electrical conduction. The partial differential equation for 
temperature distribution in a Bolid conducting heat can be compared to 
the partial differential equation for the distribution of electrical 
potential produced by a current flowing through an electrical network 
with lumped parameters. The rate at which thermal energy is stored can 
be compared to the rate of flow of charge into a capacitor. Therefore 
a resistance-capacitance electrical circuit serves as an analogue for 
heat flow in a slab. In this analogy, electrical current corresponds 
to heat flow, electrical voltage differences to temperature differences, 
electrical resistance corresponds to thermal resistance, and electrical 
capacity corresponds to thermal capacity. The physical system is lumped 
or sub-divided into geometrical sections of simple shape and form. The 
thermal properties of each section are considered to be concentrated at 
the central nodal point of each section or volume, and heat 1 b imagined 
to be conducted between nodal points through a network of fictitious 
heat-conducting rods of appropriate thermal conductance. Details of 
procedures used in the development of this analysis are given in the 
Appendix. 


3.3 CARBON MONOXIDE MEASUREMENTS 

Air samples were drawn continuously by a vacuum pump from the 
front, center, and rear of the shelter during the human occupancy tests 
(Run # 3 ) and put through the carbon monoxide measuring equipment. Car¬ 
bon monoxide samples were also drawn from these locations with the 
shelter empty while operating the gasoline driven motor generator to 
study its effect on the shelter environment. 
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CHAPTER 4 


RESUITS AND DISCUSSION 


4.1 TEMFERATJRE STUDIES 

Figure 9 shews the Bhelter temperature during Run No. 1 (simulated 
occupancy). The relationship of center shelter temperature to the air 
temperature at the entrance to the shelter and the air temperature where 
It enters the shelter proper, l.e. over the inner door, are shewn. The 
effective temperature in the center of the shelter is also shown. Ef¬ 
fective Temperature as described by ASHRAE is an empirically determined 
index of the degree of warmth perceived on exposure to different combi¬ 
nations of temperature, humidity, and air movement. The effective tem¬ 
perature index cannot be measured directly, but is determined from dry 
and wet-bulb temperatures and air motion observations by reference to 
an Effective Temperature Chart prepared by the ASHRAE. This was pre¬ 
pared from studies in which trained subjects compared the relative 
effects of various air atmospheres by occupying, alternately, two ad¬ 
joining air conditioned rooms. The acceptable effective temperature 
range for normally clothed men and women, at rest, during winter and 
summer months is 63 to 75 . In an emergency situation this upper limit 
would undoubtedly be raised. At present an effective temperature of 
85°F is considered acceptable for survival shelters. 

At the beginning of the simulated Bhelter occupancy test (Run No. l) 
the air entering the shelter was at 75 degrees and the center room tem¬ 
perature was 5 degrees lower. Within 4 hours the shelter air temperature 
was at 80° and never dropped below 75° for the duration of the run. The 
shelter temperature reflected the diurnal variations but was considerably 
damped. After the first day of occupancy the entrance air varied from 
94 to 50 degrees and the shelter temperature varied from 90 to 75 degrees. 
The effective temperature varied from 72 to 82 degrees in the center of 
the shelter which is somewhat higher than the optimum effective tempera¬ 
ture. 


The shelter temperature during Run No. 3 (Human Occupancy) is shown 
in Fig. 10. The relationship of the center shelter temperature to the 
entrance and inlet air temperature are shown. The effective temperature 
in the center of the shelter is also shown. It will be observed that 
the ambient air temperature at the start of this test was approximately 
10 degrees lower than in the simulated occupancy min. The entrance air 
reached a low of 38 and a maximum of only 67 . The minimum shelter 
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temperature was 70 and a maximum of 82. Although the temperature record¬ 
ings showed a low of 70 i n the center of the room the air entering the 
space dropped to as low as 45 degrees. Due to poor inlet air distribu¬ 
tion there were cool drafts in the front portion of the shelter. To 
minimize the air intake problem the blower was turned off for short 
periods during the coolest part of the night. This procedure was fol¬ 
lowed for the first 6 days. During the sixth day the occupants felt 
that they had corrected the difficulties by bending the baffles in the 
inlet duct; however, it will be observed on Fig. 10 that after the 
sixth day there was a coincidental rise in the inlet air temperature. 

The low entrance air temperature for the first 6 days was 36 degrees 
and on the 7th and 8th day it was 48 and 42 respectively. After the 
first day of occupancy the center room temperature fluctuated from 
70° to 82° and the effective temperature varied from Q\-° to 74°. 

Although Run No. 2 was too short to give any significant tempera¬ 
ture gradients through the soil the importance of the ventilating air 
in removing moisture from the shelter was demonstrated immediately after 
the shelter was closed down. Within l-l/2 hrs after shutting off the 
air supply, moisture was condensing on the ceiling and within 2 hrs the 
entire metal surface was covered with moisture which started running 
down and collecting at the bottom. Within 10 to 15 min after turning 
on the blower at the conclusion of the run, it was comfortable near the 
front of the shelter and within 30 to 40 min it was comfortable through¬ 
out the shelter. 


Analysis of Heat Losses in Shelter 


The heat put into the shelter by either the simulators or human 
occupants must be expelled through the shelter walls or lost to the 
ventilating air. Therefore the total heat input to the shelter equals 
the heat lost to the air plus the heat lost to the soil, or 


Q total = 


air 


+ a 


soil 


( 1 ) 


In the simulated occupancy test the heat picked up by the air consisted 
of three components: qp - the heat expended in vaporizing the added 
water and heating this vapor to the exhaust temperature; qg - the heat 
expended in raising the water vapor in the inlet air to the exhaust 
temperature; and q3 - the heat used in raising the inlet air to the 
exhaust temperature. 


Q air = *1 + *2 + * 3 


( 2 ) 
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TEMPERATURE (*F) 



1200 2400 2400 2400 2400 2400 2400 2400 2400 2400 2400 2400 2400 2400 

12/4 12/5 12/6 12/7 12/8 12/8 12/10 12/11 12/12 12/1} 12/14 12/IB 12/1 

TIME IN HOURS 


Fig. 10 Temperature Records of l4-Day Human Occupancy of Shelter. Run No 











l4-Day Human Occupancy of Shelter. 


Run Ho. 3 


Shelter human occupancy studies. 
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a j = hg at tdg - Kf at 66°P x lbs of 1^0 added 
qg = hg at tdg - hg at td^ x x lbs alr/hr 
qj = 0.24 tdg - 0.24 td 3 x lbs dry air/hr 

where hg is Enthalpy of Saturated Vapor 
tag is Exhaust Air Temperature 

hf is Enthalpy of Saturated Liquid Supplied to Simulators 
tdp is Inlet Air Temperature 

S^ is Specific Humidity of Exhaust Air (lbs of moisture/lb of 
dry air) 

0.24 tdg is the Specific Heat of the Exhaust Air lbs (dry air/hr «= 
7730 lbs/hr at 1600 cfm) 

In the human occupancy test qj_ = 0, since there 1 b no water added 
to the shelter by artificial means. 

The heat lost to the soil is expressed as follows: 


0 = k A — 

^soil A dx 


(3) 


where K is the Thermal Conductivity of the soil 

A is the Area of the Shelter exterior surface 

^ is the temperature gradient with distance from the shelter. 

With the intake air to the shelter shut off during the calibration 
runs Numbers 4 and 5, the total heat input is equal to the heat lost to 
the soil or 


Sjotal ^soll 


= h A — 
K A dx 


00 


With a known quantity of heat being put into the shelter and the tem¬ 
perature gradient in the soil determined by measurement, the value of 
k can be calculated. These calibration runs (Nos. 4 and 5) were con¬ 
tinued until the heat was being lost to the soil at a constant rate as 
determined by plotting dt vs dx. The values for dt/dx at each station 
after equilibrium was established are given in Tables 2 and 3 for Runs 
4 and 5 respectively. Substitution of values of dt/dx into equation 4 
gave a value of 1.00 for k on Run No. 4 and 1.08 on Run No. 5 or an 
average value of 1.04. 
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The heat lost to the air and to the soil was calculated for each 
2-hr interval from the above equations for both the simulated and human 
occupancy tests. The results are given in Tables 4 and 5. The averaged 
results for both runs are shown in Table 6. 

It will be noted that the heat lost to the air on the simulated 
run is 68.7 # of the total as compared to 69.2 $ on the human occupancy 
run. The calculated total heat input to the shelter on the simulated 
run is only 2.7 $ lower than the measured input. 

The total average heat loss during the human occupancy run was 
48,198 Btu per hour or 482 Btu per person. This result is compared 
with typical published values in Table 7. 

It is difficult to extend the data to other localities under vari¬ 
ous climatic conditions because of the large number of variables in¬ 
volved such as, soil conditions which affects the heat loss through 
the walls, diurnal temperature and humidity variations which affect the 
maximum, minimum, and effective temperature (E.T.), and air flow rate 
which effects the rate of heat loss from the shelter and Influences the 
E.T. The data shows that the shelter was comfortable for the entire 
l4 days of human occupancy (Run. No. 3) during which the ambient temp¬ 
erature varied between 36 and 67°F. 

The maximum E.T., however during the simulated occupancy test (Run 
No. l) was 83 with an ambient temperature variation between 42° and 95° F* 
This E.T.is a few degrees higher than the optimum maximum of 75. The 
simulated occupancy test therefore, indicates that trouble may be en¬ 
countered during the summer, in areas where the ambient temperature at 
night does not drop as low as it did during this test. 


4.2 PREDICTING THERMAL BEHAVIOR OP SHEI/EER ENVIRONMENT 

The thermal-electrical analogy is a useful tool in the prediction 
of the thermal behavior of complex heat transfer systems, but according 
to Oklahoma State University, little infoimation was available in the 
literature on its use in predicting the thermal behavior of buildings. 

The thermal electrical analogy is based upon the identity of the 
fundamental equations of heat flow within a rigid body and that 
of an electrical charge in a non-inductive circuit, in the analysis 
made at Oklahoma State University on their thermal analyzer the NRDL 
shelter interior was divided into six air nodes or cross-sections meas¬ 
uring 8 ft in depth. Node or section No. 1 was at the front of the 
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•TABLE 5 (Cont’d) 
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TABLE 7 

Metabolic Rates* of Men at Different Activities 


Activity 

Average Metabolic 
Rate 

BTU per Hour 

Sleeping 


255 

Awake, quiet 


300 

Seated, at rest 


380 

Standing, at ease 


430 

Normal Shelter Living 


482 

(Determined In these studies) 


Walking, 2 mph 


760 

Walking, 4 mph 


1400 


*A11 rates except Normal Shelter Living are frcm 


ASHRAE Guide and Data Book. 
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shelter and No. 6 at the rear. A complete report of this analyses 1 b 
presented in the Appendix. Temperature variations of Air Node No. 3 
as determined by the thermal analyzer is shewn in Pig. 11 compared to 
the measured temperature in the center of the shelter during the human 
occupancy run. It vill be noted that the temperature variations as 
determined by the thermal analyzer are 10 to 15 $> lower than the measured 
temperatures. This was a preliminary study and it is believed that 
additional work would produce better results. The analysis procedure 
offers the possibility of an inexpensive and rapid method of determining 
underground shelter environment for a wide variety of climatic conditions. 


4.3 CARBON MONOXIDE MEASUREMENTS 

During the human occupancy tests the carbon monoxide concentration 
was less than 10 ppm from 0100 to 1000 each day. After 1000 each day, 
the concentration increased over a period of 1 to 2 hours to 60 to 80 
ppm where it remained until 1600 to 1800. Thereafter, the concentration 
decreased gradually to the daily lew at about 0100. The increase- in 
carbon monoxide concentration can be explained by the increased amount 
of smoking during the daytime hours. The Navy tolerance for continued 
breathing of carbon monoxide is 150 ppm. At no time during the occupancy 
test was this level closely approached. 

Separate tests were conducted with the motor generator running to 
determine if it Increased the carbon monoxide concentration in the shel¬ 
ter.- This motor generator, which is located in the entrance-way oppo¬ 
site the ventilation fan intake is intended to provide emergency elec¬ 
trical power for ventilation and lights. The concentration of carbon 
monoxide in the shelter was considerably below the 150 ppm tolerance, 
while running the motor generator with the blower on and both intake 
and exhaust vents open. Operating the motor generator however with the 
ventilating air intakes closed; the carbon monoxide concentration rose 
considerably above this tolerance. After approximately 1-1/2 hrs it 
appeared to be leveling off at 380 bo 400 ppm as shown in Pig. 12. The 
carbon monoxide concentration dropped to a safe value within a few min¬ 
utes after opening both the intake and exhaust air vents and turning on 
the blcwer. Although the motor generator is not intended to be operated 
under these conditions this test indicates the hazard that could be en¬ 
countered if the shelter occupants had not received proper instructions. 
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Fig. 11 Measured Temperature of Shelter Compared to Temperature by Thermo-Electrical Analogy for Run Ho. 





Fig- 12 Effect of Motor Generator on Shelter Carbon Monoxide Concentration 
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CHAPTER 5 


CONCLUSIONS 


The NRDL shelter design vas proven satisfactory under the conditions 
of these tests. The temperature of the shelter remained comfortable 
during the entire l4 days of occupancy, however it is difficult to pre¬ 
dict the shelter environment for other climatic conditions. 

The shelter air temperature reflects the diurnal temperature varia¬ 
tions of the entrance air but is considerably damped. After the first 
day of simulated occupancy the entrance air varied from 50 to 9^ degrees 
and the shelter temperature varied from 7*+ to 90 degrees with an effec¬ 
tive temperature variation of 74 to 80 degrees. During the human occu¬ 
pancy test, due to poor inlet air distribution, there were cold drafts 
in the front of the shelter when the ambient temperature dropped to the 
minimum temperature of 36 ° to 42°F. The maximum entrance air tempera¬ 
ture during this test vas 67°F. The air in the center of the shelter 
varied from 70 ° to 82° and the effective temperature varied from 64 to 
74 degrees. 

The major portion of the heat generated in the shelter was removed 
by the ventilating air. In both the human occupancy test and the test 
run using human simulators, 69 $ of the heat vas carried off by the 
ventilating air and only 31 $ vas lost to the soil surrounding the 
shelter. 

The 1600 cfm of air flow through this shelter reduces the carbon 
monoxide from smoking to a concentration far below the Navy tolerance 
for submarines of 150 ppm, and it exceeded the maximum allowable concen¬ 
tration (MAC) established for industry of 100 ppm for only approximately 
1 hour. The industry MAC of 100 ppm was established as maximum for 8 
hours exposure. Prom 0100 to 1000 each day the concentration averaged 
10 ppm. The peak of 60 to 80 ppm was reached each day by 1200 where it 
remained until it began to drop off between l 600 to 1800 . 

The location of the motor generator will result in dangerous quanti¬ 
ties of carbon monoxide being created if this generator is operated by 
mistake with the ventilating fan on and the intake air vents closed. 

Thermal-Electrical analyzer appears to offer practical method of 
predicting shelter environment. 



CHAPTER 6 


RECOMMENDATIONS 


It is recommended that! 

(1) The position of the motor generator be changed or a safety 
measure be devised to prevent it from being operated at conditions 
other than those for which it is intended. 

(2) Changes be made in the inlet ducting to eliminate cold drafts 
at the center of the shelter during cold weather. 

(3) Consideration be given to the further development of Thermal- 
Electrical Analogy as a tool to predict shelter environment over a 
wider range of climatic conditions. 
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Subscripts 


- Capacity (Thermal and Electrical) 

- Specific Heat ® Constant Pressure, BiO/#°F 

- Voltage. Volts 

- Current - Amperes 

- Thermal Conductivity. Btu/hr.ft.°F 

- Mass Flow Rate, fftn/hr 

- Quantity of Heat Energy per Unit Time. Btu/hr 

- Resistance (Thermal and Electrical) 

- Density, ^m/ft^ 

- Temperature. °F 

- Time 

- Volume, ft-* 

- Capacitance (Electrical & Thermal) 

- Diode 

- Ground Connection 

- Resistance (Electrical & Thermal) 

- Electrical Circuit 

- Thermal Circuit 
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CHAPTER I 


INTRODUCTION 


The passive R-C electrical analogy for thermal systems 
is a powerful tool in the prediction of the thermal behavior 
of complex heat transfer systems, and many investigations 
using this method oan be found in the literature. However, 
very few applications to the study of the thermal behavior 
of buildings can be found. 

This paper presents a thermal analysis using the 
electrical analogy to determine the temperature variation 
within a one hundred man underground shelter. The system 
studied is a Navy Ammunition hut of steel-arch type structure, 
48 feet long, with an arch radius of 12.5 feet, which is 
covered by a minimum of three feet of earth. 

This structure was built by the United States Naval 
Radiological Defense Laboratories, (USNRDL), approximately 
forty miles from San Francisco, California, where physical 
and analytical studies of this system are in progress. 

These studies are directed toward determination of environ¬ 
mental conditions within an underground shelter, occupied 
by one hundred men for extended periods of time. Rough 
drafts of these studies we^e presented to the Mechanical 
Engineering Department of the Oklahoma State University by a 
representative of USNRDL, In the analysis presented, the 
system and boundary conditions are the same, to the extent 
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that it is possible, as in the USNRDL tests. This permits 
the results of this analysis to be compared with those of 
the USNRDL tests. 

In an analysis such as this, the main problems are not 
associated with the thermal-electrical analogy nor with the 
thermal analyzer. It is one of representation. The actual 
physical conditions and boundary conditions must be represented 
by an equivalent thermal network or circuit. In accomplishing 
this, the physical system is lumped or subdivided into a number 
of geometrical subvolumes. The thermal properties of each 
subvolume is considered to be concentrated at the central 
nodal point of each subvolume, and heat is imagined to be 
conducted between nodal points through a network of fictitious 
heat-conducting rods of appropriate thermal conductance. It 
should be apparent that the lumping of the physical system 
and the choice of values of properties assigned to nodal 
points will greatly influence the accuracy of such an analysis. 
Methods must also be devised for the representation of boundary 
conditions and any other parameters which might affect the 
thermal behavior of the system. Chapter II presents the 
equivalent thermal network that represents the system and this 
is related to an electrical network in Chapter III by means 
of the thermal-electrical analogy. The electrical circuit is 
then programmed on the thermal analyzer which leads to an 
analysis which is interpreted in terms of thermal properties. 
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CHAPTER II 


THERMAL CIRCUIT 

In an analysis such as this, it is necessary to idealize 
or simplify the problem so that it may be represented by a 
manageable thermal circuit. In doing this the following 
assumptions were made. 

(1) The shelter is considered to be made up of two ends, 
floor, and siding, all of which represent paths of 
heat transfer from the enclosure. 

(2) All thermal energy transferred through the structure 
is unidirectional. This is to say that there i8 

no heat transferred between the paths mentioned 
in (1). 

(3) Physical conduction pathB are replaced by sections 
or lumps which represent the thermal properties of 
the path. 

(4) The thermal properties of the lumps are constant 
and are mean values over the temperature range 
encountered. 

(5) Net radiation exchange within the enclosure is 
negligible, except that which is included in the 
conductance from the men to the walls. 

(6) Solar radiation input into the system is 
considered negligible. 
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(7) Phase changes, such as condensation on the walls 
are not considered. 

(8) Contact resistances through the conduction paths 
are negligible. 

The interior of the enclosure is divided into six nodes 
as pictured in Figure 1. The paths of heat transfer from 
these nodes are shown in Figures 14, 15 and 16, in the Appendix. 



Figure 1. Overall Schematic of the Enclosure 
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The walls, adjacent rock, and earth are assumed to form 
a series resistance to the radial exchange of energy between 
the interior nodes of the shelter and its surroundings. The 
floor and adjacent earth also represent series resistances to 
the energy exchanges. In addition, the end walls of the 
enclosure represent additional conduction paths to be added 
to each end node. 

An interior node can be represented by the network as 
shown in Figure 2. The resistances represent thermal resistance, 
and the capacitors represent thermal capacity of each lump in 
the conduction path.,, Figures 14, 15 and 16, will aid in the 
understanding of the physical picture. 



Conduction Path Between 
Internal Node, Floor and Earth 


Internal Node 


Conduction Path Between 
Internal Node, Siding and Earth 


Figure 2. Thermal Circuit Representation of an Interior Node 



The boundary conditions imposed on the thermal system 
consist of a heat input into the shelter to be maintained 
at 470 B/hr-man for one hundred men, and 1600 CFM of 
ventilating air forced in at one end of the shelter. 

The one hundred men can be represented in the thermal 
circuit by a heat generation or energy input, "q* 1 into each 
internal node of the enclosure. 

The ventilating air presents a problem in that there is 
no analogy such as exists for the conduction of heat. 
Physically the problem is as followst 
The ventilating air is assumed to vary between 40°F 
and 60°F, while the enclosure temperature is expected to be 
ten to fifteen degrees above the mean temperature of the 
ventilating air. The flow of cooler air through the 
enclosure establishes a temperature gradient with the 
temperature increasing in the direction of the air flow. 

This gradient affects a flow of energy in a direction 
opposite to that of the air flow. This is shown schemat¬ 
ically in Figure 3. Another important consideration is 
that the temperature of a particular node is affected by 
conditions upstream of the node only. 

_Air Flow_w 

O O O O O 0 o 

(0) (1) (2) (3) (4) (5) (6) 

_ Heat Flow 

Figure ). Schematic Showing Direction of Heat and Air Flow 

in the Enclosure 
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The incominp, air impresses its temperature or potential 
on its entrance point which will be designated node (0). The 
increase of energy of the air moving from (0) to (l) can be 
described by q^ - ^ a assistance is defined as 

R ~ 1/mCp, then Figure 4 would be an analogous thermal circuit* 
The 8ymbol-| | | represents the potential of the incoming air. 



Figure 4. Ventilating Air Analogy Applied to the First Air Node 

In considering node (2) which is only dependent on node 
(1), while node (1) is not dependent on node (2), we haves 

^2 = mc p (T2 - Ti) 
or 

q 2 = mc p (T 2 - t 0 ) - m C p (T x - t Q ) 
and 

32 = mc p (^2 " ^0^ ~ 3]_ 

This is represented schematically in Figure 5. 
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Figure 5. Ventilating Air Analogy Applied to First and 

Second Air Nodes 

The symbol—— in Figure 5 is an electrical symbol 
representing a diode. The diode used a9 shown, prevents 
the flow of energy from (l 1 ) to (1). This is a necessary 
condition in that node (l) is independent of any point 
down stream of it. Because the resistances in Figure 5 are 
equal, the potential of node (l 1 ) will equal that of node (l), 
It follows that q^ dependent only upon R» tg and t^. 

The circuit is now evaluated to determine if it is 
analogous to the physical conditions it represents. 

Physically the exchange of energy at node (2) is: 

q 2 out = mCpt 2 

q 2 in - mCpt^ 
and 

q 2 net = mCp(t 2 - t^) 

It is concluded that the physical conditions at node (2) 
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are satisfied by the thermal circuit of Figure 5. It is 
pointed out again that the use of the diode prevents the 
flow of energy from (1') to (1), An energy balance on 
node (1) in Figure 5 would be identical to an energy 
balance made on node (l) in Figure 4. 

The circuit for all six nodes is shown in Figure 6. 



Now evaluating the enclosure physically? 

Energy into enclosure from air flow q^ n = mc p tQ 
Energy out of enclosure from air flow q ou t m mCpt^ 

The net flow out is q net “ mc p (t6-to) 

Summing the flow of energy from the enclosure using the 
thermal circuits 

% 8 mc p (t 6 - t 5 ) 
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q<j = mCp (t 5 - t 4 ) 
q 4 = mc p (t 4 - t 3 ) 
q 3 * mc p (t 3 - t 2 ) 

q 2 " mc p ^2 ‘ ^ 

q l s mc P ^1 “ fc o) 

q = mcp (t^ - t 0 ) 

Therefore this circuitry is truly analogous to the 
physical conditions as established. 

Figure 7 shows the thermal circuit for an interior node 
with boundary conditions included. 
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Figure 7. Thermal Circuit Representation of an Interior Node 
with Boundary Conditions Included 


CHAPTER III 


ANALAQOUS THERMAL CIRCUIT 

The thermal electric analogy! which is based upon the 
identity of the fundamental equations of the flow of heat 
within a rigid body and that of a charge in a non inductive 
circuit, is well known.^ Figure 8 is therefore a circuit 
diagram, representing the electrical as well as the thermal 
cirouit. The scale factors or ratios derived to complete 
the analogy are shown in Table I. 


TABLE I 

SCALE FACTORS USED IN THE ANALOGY 


Quantity 

Units 

Scale Factors 


Thermal: 

Electrical: 

msssm 

Value 

Time 

Hrs. 

Sec. 

? IT, 

1 hr/sec 

Capacity 

Btu/°F 

Farads 

C 7 lc e 

10 8 B/°F 
farad 

Resistance 

°F/Btu/hr 

Ohms 

/? 7 /*e 

10“ 8 hr°F/B 

Potential 

Op 

Volts 

tie 

1 °F/volt 

Rate of Energy 
Transfer 

Btu/hr 

Coulombs/sec 

or 

Amperes 

9U 

10 8 Btu/hr 
amp 


^-The writer has written a paper on this subject. 
Bibliography Reference No. 8. 
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Of the first three ratios in Table I, thermal-electrical 
analogy permits the selection of two. The third one being 
fixed by the selection of the other two. The relating 
equation is 7?/% = ^/C e * ^//Qe 

It was desirable in this analysis to choose a time ratio 
wherein days would be sealed down to a matter of minutes. 

For convenience, one hour of thermal time was set equal to 
one seoond of eleotrical time. In ohoosing the capacity and 
resistance scale factors, the limitations were the magnitude 
of the electrical resistors and capacitances available. It 
was also desirable to keep capacitor leakage at a minimum. 

The choice of these scale factors was therefore a compromise. 

Of the last two ratios in Table I, the selection of one 
would determine the other. The relating equation for this is; 

T /f = 9// X*T//? e 

Since R t /l 0 had been set * choice of T /E and qA Was 
again a compromise between selecting values of current of an 
order of magnitude greater than that of the anticipated 
leakage current and voltages that were readily available. 

The boundary conditions imposed upon the thermal circuit 
must be duplicated with regard to the electrical circuit. As 
mentioned previously, the diurnal temperature variation of the 
ventilating air was assumed to be between 40°F and 60°F. 

Having selected T/E equal to l°F/volt, this was represented 
by a Donner Analog Computer, programmed to produce an output 
voltage which varied sinusoidally between 40 and 60 volts, 
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with a 24 second period. According to the time ratio, 24 
hours is represented by 24 seconds. 

The input from the one hundred men is 47,000 B/hr. 

For each air node this amounts to 7,833 B/hr, or in the 
electrical circuit 78.3 yU amps. This was accomplished 
with a 400 Volt D.C. power supply connected by a large 
resistance (approx. 4 meg._/--\_ ) to each air node. These 
resistances were adjusted while the problem was being run, 
so that proper current flow was achieved. The resistances 
labeled 2R, 3 R, 4R and 5R, were also adjusted during running 
to insure that the desired voltage duplication was achieved. 

Figure 9 and 10 are photographs of the test setup. The 
problem board upon which the electrical circuit is arranged, 
is the panel on the right side in Figure 10. The capacitor 
leads are seen on the left in the same photograph. 
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CHAPTER IV 


MEASUREMENTS 


During the course of this study, voltage measurements 
were taken at every node in the circuit. These measurements 
were recorded on a Sargent Recorder which received its input 
from the output of a very high impedence voltage measuring 
instrument or electrometer. The electrometers are an integral 
component of the Mechanical Engineering Department thermal 
analyzer. This data was not presented in this report because 
it was taken only as a means to determine if the circuit was 
performing properly. 

Data presented in this report is in the fom of graphs 
and is the temperature time variation of the two middle air 
nodes. The lower left hand portion of the curves, (Figures 
11 and 12), are the readings at the middle nodes prior to 
the entrance of the hundred men into the shelter. What 
appears to be a discontinuity in the curve is the application 
of the input representing the one hundred men. This was 
included so that the temperature buildup with the entrance 
of the men could be shown. Values for the lower left hand 
portion of the curve should not be considered valid, since 
the ventilating air analogy was based upon the assumption 
that the ventilating air was of lower temperature than the 
air inside the shelter. 

Measurements of the leakage resistance of each node in 
the circuit was recorded. Current leakage in the thermal 
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analyzer represents an effective loss of energy from the 
thermal system and is therefore a source of error. Although 
this is undesirable, it must be realized that it exists. 
Measurement of the leakage resistance was accomplished by 
charging the capacitance on each node to 80 volts and 
observing the time for the node to drop to a given potential. 


56 



9 ( 



Figure 11. Temperature Variation of Air Node No. A Compared 

with USNRDL Data 
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Figure 12. Temperature Variation of Air Node No. 3 Compared 

with USNRDL Data 
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Figure 13. Comparison of Incoming Air Temperature 
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CHAPTER V 


RESULTS AND CONCLUSIONS 

The results of this study are the curves shown in 
Figure 11 and Figure 12, Also shown on these figures is a 
plot of the data representing the center room temperature 
during the USNRDL tests. Comparison of the center room data 
can best be made when our boundary conditions make the best 
possible fit with actual conditions as experienced by USNRDL. 
Figure 13 is a plot of the actual temperature variation of 
the ventilating air during the USNRDL test compared with ^ 
our assumed temperature variation. 

The difference in the results of this study and the 
USNRDL test is of the order of to 1C® and is very good 
when one is considering a complicated heat transfer system. 

When a differential equation is replaced by a finite difference 
equation, errors result. These errors are difficult to 
determine, but for the data taken in this study they are 
estimated to be negligible, especially with respect to 
leakage errors. 

For the two nodes of interest, namely the two middle 
air nodes, the average leakage resistance in the conduction 
paths from these nodes is approximately 230 meg ohms. This 
corresponds to a leakage current of approximately 1.5 /U amps. 
This is very small compared to the 78 yO amps into each of 
these nodes. It is also small compared to the current which 
is flowing out of the nodes representing the ventilating air, 
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which is of the order of magnitude of 100 p, amps. 

A more accurate comparison with the USNRDL tests is not 
possible because the writer does not have detailed information 
as to the instrumentation and techniques used in temperature 
measurements during the USNRDL tests. The location of 
thermometers and thermocouples with relation to heat sources 
and proper shielding is very important. Any deviation from 
the proper technique could easily account for the difference 
shown. # 

In a last analysis, it may be stated that a thermal 
circuit, solution can only be as good as the representation 
of the actual physical system and the ability to evaluate 
the circuit parameters and boundary conditions. This 
study has demonstrated the value of the thermal analyzer in 
an economical and rapid solution of a complex heat transfer 
system. 


^Author's note: Mr. Othling did not have complete experi¬ 
mental details at the time this report was prepared. 
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APPENDIX. * 


COMPUTATION OP THE PARAMETERS 


Thermal Resistances and Capacities in the Conduction Path 
Through the Siding (Figure 14) 

Material 8 gage sheet steel 

Thickness .1644 inches 


e 

k 


488 #/ft 3 
27 Btu/hr.ft«°P 
.112 Btu/#°F 


Properties of the earth adjacent to the shelter! 


f 

k 


108 #/ft 3 

.705 Btu/hr.ft„°F 


C p .26 Btu/#°F 


The earth properties are those used in the USNRDL 
analytical study and are weighted values determined 
according to procedures used by the Army Corps of 
Engineers. 

Properties of the Air Space: 
t 13.7 ft 3 /# 

C p ,24 Btu/#°F 


■^Appendix to Othling's original document. 




For reflective surfaces with an effective emissivity of 
0.05, the "Heating,, Ventilating and Air Conditioning Guide" 
recommends an inside surface film coefficient of 0.8. 

This includes convection and radiation effects and is 
independent of the position of the surface and the 
temperature difference. Even though this particular 
choice of "h" is subject to discussion, it is of minor 
importance, in that it represents a small percentage of 
the total resistance to heat flow through the walls of 
the shelter. 


6k 


h - .8 Btu/hr.ft.°F is the value used in the USNRDL 

tests and likewise will be used in this study, 

n = _Z_ r _'_ 

K 1 fi/j .Q*B'X/*.S*6 

- 3.98x10“^ hr°F/Btu 
r r jqCoY - .-jL?.* %JL '*<* . *,# 

c i r “ /3.7xz 


= 34.4 Btu/°F 


R ~ /G/A) f nr . atna n r* /■» 

*2 2 PA/. ror 8ma11 r 2 /r l . 

- s/ZVjLste = .eowoArOF/Btv, 

„ - AC*/ a #46 X.//*X./**¥XPx/UA* 

^2 “ f ' 3 /£ 


= 236 Btu/°F 


*3 


2rrr/. . 7<>j'x fT/a 


= .00434 hr°F/Btu 


C 


3 



/#£ *.2*xF/**-/:*)<9 

~~S~ 


- 1896 Btu/°F 


Thernial Resistances and Capacities in the Conduction Path 
Through the Ends of the Shelter (Figure 15) 

Material 7 gage sheet steel 

Thickness .1793 inches 
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Figure 15. Conduction Path through End Walls 



66 





k .75 Btu/hr.ft. C< P 

C p .2 Btu/#°F 



)Air Node 

JR? 

► 

11 

:r 8 

3 Floor 

7 

< 

Si 



■ 

ni-h^ 
_ £ 

1*9 Earth 

3 

>' - 

; 


- *r* 




mmgm 

i 1 


Figure 16. Conduction Path through Floor 


p - -jL. 
/M 

*8 ■- Si- 

Cg = Pfr? 

"»> = itfr 

c, = /Of*' 


= 1 = .00625 hr°F/Btu 

78x25x8 

= 1.5 = .833xlO~ 3 hr°F/Btu 

12 x.75x25x8 

= 400x.2x25x8x3 = 4,000 Btu/°F 
12 


= 1 = .00708 hr°F/Btu 

7705x25x8 

=108x.26x2x25x8 = 11,200 Btu/°F 


Resistance R = l/mCp defined by ventilating air analogy; 
Flow rate = 1600 C.F.M. 

R * l/mc = 13.7 = .594xl0" 3 hr°F/Btu 

p l600x.24x60 
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TABLE II 


SUMMARY OF RESISTANCES 


Symbol 

R t 


\*h 

399xl0“ 6 

39„9/’^~\_ Convective & Side Walls 

h'h 

.00434 

.434 meg^-s. Side Walls & Earth 


.0509 

5.09 meg./~\-Convective & End Walls 

H 5 f % 

.00578 

.578 End Wall9 & Earth 

B, t % 

.00708 

.708 meg.'"-'^ Convective & Floor 

Rg / R 9 

.00791 

.791 meg./^- Floor & Earth 

R = l/mc p 

594xlO“3 

59.4 /(ss.. Ventilating Air 


X R = Rt * 

10" 8 hr°F/B 



ohm 
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TABLE III 


SUMMARY OF CAPACITANCES 


Symbol 

5 

c e 


Cl 

34.4 

.34 l/U/ 

Internal Air 

°2 

236 

2.3^/ 

Side Wall 

°3 

1896 

18.9^/ 

Earth Adj. to 
Side Wall 

C 5 

200 

2yL// 

End Walls 

c 6 

1.38x10^ 

13^/Z 

Earth Adj. to 
End Wall 

C 8 

4,000 

UO/// 

Floor 

c 9 

11,200 

112/// 

Earth Adj. to 
Floor 


*c 8 c t " 

10 8 B/°F 

Farad 
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